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Abstract: We present a novel technique for three-dimensional (3D) 
image processing of complex fields. It consists in inverting the coherent 
image formation by filtering the complex spectrum with a realistic 3D 
coherent transfer function (CTF) of a high-NA digital holographic micro- 
scope. By combining scattering theory and signal processing, the method 
is demonstrated to yield the reconstruction of a scattering object field. 
Experimental reconstructions in phase and amplitude are presented under 
non-design imaging conditions. The suggested technique is best suited 
for an implementation in high-resolution diffraction tomography based on 
sample or illumination rotation. 
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1. Introduction 

High-resolution three-dimensional (3D) reconstruction of weakly scattering objects is of great 
interest for biomedical research. Diffraction tomography has been demonstrated to yield 3D re- 
fractive index (RI) distributions of biological samples [1]. For the use of such techniques in the 
field of virology and cancerology, a spatial resolution in the sub-200nm range is required. Con- 
sequently, experimental setups must shift to shorter wavelengths, higher numerical apertures 
(NA) and steeper illumination and/or sample rotation angles. However, the scaling of resolu- 
tion to high-NA systems introduces strong diffraction and aberration sensitivity [2]. The use 
of microscope objectives (MO) under non-design conditions, in particular for sample rotation, 
introduces additional experimental aberrations that may further degrade resolution [3,4]. Un- 
fortunately, the theory of diffraction tomography cannot correct for these conditions since it is 
based on direct filtering by an ideal Ewald sphere [5]. 

We present a novel technique that reconstructs the object scattered field using high-NA MO 
under non-design imaging conditions. Opposed to classical reconstruction methods like filtered 
back projection [6], our approach is based on inverse filtering by a realistic coherent transfer 
function (CTF), namely 3D complex deconvolution. We expect this technique to lead to aber- 
ration correction and improved resolution [7, 8]. By combining the theory of coherent image 
formation [9] and diffraction [5], the deduced theory enables reconstruction of object scattered 
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field by inverse 3D CTF filtering. Under sample rotation, experimental evidence of this theory 
is presented. It confirms that the realistic 3D CTF can be directly reconstructed from a hologram 
acquisition of a complex point source [10, 11]. For this purpose, DHM's feature of digital re- 
focusing is exploited [12]. As simulations have shown, digital refocusing bears the capacity of 
optical sectioning [13]. By regularizing the complex filter function, phase degradation by noise 
amplification is suppressed as anticipated from intensity deconvolution [14]. The effectiveness 
and importance of the proposed method is demonstrated with experimental applications. 

Independent of the non-design imaging condition, our method serves to reconstruct the scat- 
tered complex object function from a single hologram acquisition. One reconstruction does not 
feature tomography but rather optical sectioning of one rotation angle measurement. However, 
the reconstruction can be applied to various complex field acquisitions, which can be ultimately 
used for super-resolved diffraction tomography. 

2. Description of method 

The proposed method consists of three major inventions. First, for inverse filtering the three- 
dimensional deconvolution of complex fields is formalized by complex noise filtering. Sec- 
ondly, based on single hologram reconstruction, an experimental filter function is defined. 
Third, in a rigorous approach, the filtered field is used to retrieve the scattered object func- 
tion. 

2.7. Regularized 3D deconvolution of complex fields 

For a coherently illuminated imaging system, the 3D image formation of the complex field U 
is expressed as the convolution of the complex object function, called o, and the complex point 
spread function (APSF), called h [9]: 



r2-ri)dxidjidzi, 



(1) 



where r = (x, denotes the location vector in the object space f\ and the image space r2 as 
shown in Fig. 1(a). Equation (1) can be recast into reciprocal space by a 3D Fourier transfor- 
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Fig. 1. Optical transfer of a point source in real and reciprocal space. In scheme (a), a 
practical Abbe imaging system with holographic reconstruction. In scheme (b), full Ewald's 
sphere under Born approximation in the reciprocal object plane according to Eq. (11). 

mation ^ defined as: 



JJ (r2) exp[i27r(^ • r2)]dx2d}72dz2. 



(2) 



The reciprocal space based on the free-space (index of refraction n = 1) norm of wavenumber 
k with wavelength A, relates to the spatial frequency v and wave vector I = (k^^ky^k^ by 
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2k 

k=\k\=2KV = —. (3) 
According to the convolution theorem, applying Eq. (2) to Eq. (1) results in: 

^{U{f2)} = ^{o}^{h}. (4) 

G(k) 0(k) c(k) 

Conventionally, the 3D Fourier transforms of U, o, and h are called G, the complex image 
spectrum, O, the complex object spectrum, and, c, the coherent transfer function (CTF), as 
summarized in Eq. (4). The latter is bandpass limited through h, with the maximal lateral wave 
vector, 

kxy^c = k sin a, (5) 
and the maximal longitudinal wave vector 

^z,c = ^(1 — cosa), (6) 

scaled by k from Eq. (3). The angle a indicates the half-angle of the maximum cone of light 
that can enter into the MO (cf. Fig. 1) and is given by its NA = /t/sina (rii is the immersion's 
index of refraction). Through Eq. (4), the complex image formation can be easily inverted: 

^(n) = 0(^)exp[-i27r(^-ri)]dMM^. = ^"'{^}- (7) 

The 3D inverse filtering can be directly performed by dividing the two complex fields G and 
c. As known from intensity deconvolution [15], the inverse filtering method in the complex 
domain suffers from noise amplification for small values of the denominator G(k)/c(k), partic- 
ularly at high spatial frequencies. 

As stated by Eq. (4), the recorded spectrum G(k) is physically band limited by the CTF, 
thus it can be low-pass filtered with the maximal frequency kjy^c of Eq. (5) in order to suppress 
noise [8]. However, small amplitude transmission values within the band pass of the 3D CTF 
may still amplify noise. The noise amplification results in peak transmission values in the de- 
convolved spectrum, which add very high modulations in phase. Thus, phase information could 
be degraded through amplitude noise. To reduce noise degradation effectively we propose a 
threshold in the 3D CTF of Eq. (7), such as: 



^(^) = \ ; r Ml ' I r • (8) 
I 1 •exp[i-arg[c]] if|c|<T 

For modulus of c smaller than T, the CTF's amplitude is set to unity, so that its noise ampli- 
fication is eliminated while its phase value still acts for the deconvolution. By controlling T, 
truncated inverse complex filtering (t << 1) or pure phase filtering (T=l) can be achieved. 
Therefore, the deconvolution result depends on the parameter T. Compared to standard regu- 
larization 3D intensity deconvolution [14], the threshold acts as a regularization parameter in 
amplitude domain while the complex valued domain is unaffected. The influence of this param- 
eter is discussed in section 3. 

Note that thresholding based regularization is comparable to other schemes, like gradient 
based total variation regularization [16]. However, the presented method does not assume spar- 
sity, as known from compressive sensing approaches. Instead, it is based on the inversion of the 
coherent image formation of Eq. (1). 
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2.2. 3D field reconstruction of 2D hologram 

Typically, the 3D image of a specimen is acquired from a series of 2D images by refocusing the 
MO at different planes [14]. In the proposed technique, however, the complex fields are pro- 
vided by digital holographic microscopy (DHM) in transmission configuration [17] as shown 
in Fig. 2(a). 




Fig. 2. Experimental configuration, (a) Optical setup: LD, laser diode; BS, beam split- 
ter; M, mirror; DL, delay line; SF, spatial filter; ND, neutral density filter; L, lens; TL, 
tube lens; CL, condenser lens; MO, microscope objective; RS, rotatable specimen, (b) RS 
with complex point source in MO design and non-design conditions. Insert: SEM image of 
nano-metric aperture (0 ^15nm) in aluminum film at 1 50000 x magnification, (c) RS with 
objects (see section 3) in experimental conditions with incident light along optical axis. 



Thus, the amplitude A (r2) as well as the phase 0(r2) of the hologram ^ can be reconstructed 
by convolution [12]: 



U{r2) =A(r2)-exp[iO(r2)] 



exp(iM^ J 
id'X 



/£^^(4,J2)exp [{x'2-X2f^{y'2-y2)]^ 



where r2 is a spatial coordinate in the hologram plane, and d is the hologram reconstruction 
distance as shown in Fig. 1(a). Using digital refocusing, a pseudo 3D field can be retrieved by 



varying the reconstruction distance dz2 = d^M^zi scaled by the MO's longitudinal magnifica- 
tion of = (/2//i)^. Note that opposed to MO refocusing, the physical importance of digital 
refocusing is related to the camera's distance from the system's image plane. 

2.3. Experimental pseudo 3D APSF 

The coherent imaging system can be experimentally characterized by a complex point source, 
shown in Fig. 2(b). It consists of an isolated nano-metric aperture (0 ^75nm) in a thin opaque 
coating (thickness=100nm) on a conventional coverslip [11]. The aperture was fabricated in 
the Center of MicroNano-Technology (CMI) clean room facilities by focused ion beam (FIB) 
milling in the evaporated aluminum film (thickness=100nm). For a single point object o{fi) = 
6(ri), the image field U§{f2) is the APSF h{f2). This approximation yields aperture diameters 
0 << dmin (dmin' limit of resolution), and its imaged amplitude and phase have been shown to 
be characteristic [10, 18]. 
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The coverslip is mounted on a custom diffraction tomography microscope [19] based on 
sample rotation and transmission DHM as shown in Fig. 2(a). The setup operates at A = 680^m 
and is equipped with a dry MO for long working distance of nominal NA = 0.7 and magnifi- 
cation M = 100 X . The sample rotation by 6, as depicted in Fig. 2, introduces non-design MO 
conditions of imaging [4], consecutively discussed in section 3.1. The optical path difference 
between corresponding rays in design and non-design systems can be caused by use of a cov- 
erslip with a thickness or refractive index differing from that of the design system, the use of 
a non-design immersion oil rii, the defocus of the object in a mismatched mounting medium 
[3]. Consequently, a system's defocus must be avoided through sample preparation, as dis- 
cussed in section 3. Since a dry MO is used with matched coverslips, we expect the tilt to 
introduce the main experimental aberration, apart from MO intrinsic aberrations. 

In order to demonstrate the importance of the proposed technique to diffraction tomography 
by sample rotation, experimental holograms are recorded for tilt positions as well. The digitally 
reconstructed pseudo 3D APSF are shown for two positions in Figs. 3(a)-3(b), without tilt and 



Fig. 3. Measurement of complex point source in MO design and MO non-design imaging 
conditions (cf. Fig. 2). The experimental APSF sections in (a) yield design MO imaging 
conditions (6=0°), whereas sections in (b) yield non-design conditions (0 = 15°). The left 
side images show \h\ central sections and the right images arg[/z], respectively. Colorbar, 
Scalebars: 2/im. 

In MO design conditions, the complex field in Fig. 3(a) indicates a typical point spread 
function pattern [10]. In amplitude, the diffraction pattern is similar to the Airy diffraction 
pattern, derived from the Bessel function of first kind Ji . The phase part oscillates at Ji 's roots 
with spacing A /NA from — tt to k [7]. Nonetheless, the axial sections in Fig. 3(a) are prone to a 
spherical like aberration due to intrinsic aberrations [2]. The lateral sections show a good axial 
symmetry, hence the absence of strong coma or astigma like aberrations. 

In the case of 6 = 15°, the field in Fig. 3(b) features asymmetries of the diffraction pat- 
tern in amplitude and phase. The aberration can be especially well observed in a lateral phase 
distortion. Likewise, the asymmetric aberration is also expressed in the axial direction of the 
APSF. The introduction of coma-like aberration [2] is due to the tilted coverslip system. Just 
like above mentioned non-design conditions [3], the tilt results in optical path differences in the 
experimental system, which act as an additional aberration function [20] . 

The dependence of the APSF's amplitude A§ and phase on the sample rotation (cf . Fig. 2) 
is defined by 
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for e = 15°. 




with the illumination wavevector k'^ = {kQ^x-,kQ,y,kQ^z) the laboratory reference frame, mean- 
ing relative to the optical axis. In the demonstrated case of sample rotation, k'^ = (0, 0, 1) relative 
to the optical axis does not change. However, the illumination relative to the sample does change 
and the incident field vector can be expressed as ^0 = (O,sin0,cos0) in the sample frame of 
reference, according to Fig. 2(b). 

The complex point source technique allows registering the scattered h background illumi- 
nation free. As a result, the SNR is advantageous and the required h can be directly used for 
Eq. (2). The amplitude and phase of the recorded field (r2) and (r2) then correspond 

is) ->■ (s) ->■ 

to the scattered components A\L (r2) and (r2). 

Eventually, it is the design of the complex point source (cf. Fig. 2) that guaranties a high 
accuracy of the APSE with respect to realistic imaging. The complex point source is directly 
located on a conventional coverslip, also used for object imaging. Thus, the APSE acquisition 
does not require any modifications of the imaging system, is stable, and most importantly, it is 
representative for object imaging under the same conditions. Even, for defocused objects, the 
experimental APSE can be estimated by flipping the complex point source and immersion it 
with approximated defocus. 

2.4. Experimental 3D CTF 

If the incident illumination field on the scatterer is in direction ^0 and the scattered field is 
measured in direction of k, the first-order Born approximation states that the 3D CTE is given 
by the cap of an Ewald sphere defined for 

K = k-ko, (11) 

as schematically shown in Eig. 1(b) for DHM [21]. NA limits the sphere, so that only part of the 
diffracted light can be transmitted. The experimental DHM's 3D CTE can be directly calculated 
by applying Eq. (2) to Eq. (10) and is depicted in Eig. 4. 

Erom this reconstruction, the NA can be directly measured through the subtended half-angle 
a according to the cut-off frequencies of Eqs. (5) and (6). The experimental CTE includes 
experimental conditions, i. e. aberrations, as well. By comparing Eig. 4(a) and Eig. 4(b) the 
impact of illumination becomes apparent. Due to the sample rotation, the MO can accept higher 
frequencies from one side of the aperture, while frequencies from the opposed side are cut off. 
As a result, the CTE is displaced along the Ewald's sphere [22]. Note that this displacement is a 
combination of translation and rotation if the rotational center is not coincident with the sample 
geometrical center. 

Thus, the 3D CTE can be written as a function of K: 



c{k-ko)=A^s\^-^o)-^W i^'s'{k-ko) 



(12) 



where a hat indicates the Eourier component in amplitude A and phase O as summarized in 
Eig. 5(a). Similarly to the 3D CTE reconstruction, the three-dimensional complex spectrum 
G(k — ko) is calculated from experimental configurations shown in Eig. 2(c). 

2.5. Scattered field retrieval 

In the case of transmission microscopy, the APSE in Eq. (1) is not directly convolved with the 
complex object function o. According to diffraction theory [5], the total field o can be expressed 
as the sum of the incident field o^^ in direction of k^ and the scattered field o^^'. 
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a =35° 




a =42.5° 



Fig. 4. Experimental 3D CTF in different imaging conditions. The experimental CTF in (a) 
obtained for MO design imaging conditions (6 =0°), whereas the CTF depicted in (b) is 
obtained for non-design conditions (6 = 15°), according to Fig. 2. The upper row shows 
the top view on the CTF and bottom row shows the side view through the CTF for kx = 0, 
respectively. Colorbars. 



(13) 



where 



oi:\ri)=Ai:\ri)-txpU^\ri) 

ko k) Ik) 



(14) 



(s) ->■ (s) ->■ 

with the scattered field amplitude Ai Mri) and phase Oi (n). On substituting from Eq. (13) 

ko ko 

into Eq. (1) and using Eq. (7) we see that the complex deconvolution satisfies 



c(k-kQ) 



(15) 



The subtracted convolution term in the numerator of Eq. (15) can be identified as c(k — 
ko)^{o^^ iri)}, the reference field of an empty field of view [23]. Suppose that the field in- 

ko 

cident on the scatterer is a monochromatic plane wave of constant amplitude A^^\ propagating 
in the direction specified by ^o- The time-independent part of the incident field is then given by 
the expression 



exp 



i^on 



(16) 



and Eqs. (15) and (16) yield 
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Fig. 5. Scheme of reciprocal space. In image (a), a ID coherent transfer function as given by 
the complex point source [see Eq. (12)]. In image (b), image's spectrum with background 
illumination ko ^ 0. 



^{4o(n)} = 



G(X-Xo)-A« 5(^-^0) 



(17) 



On the other hand, according to Eqs. (13), (14), and (16) the image spectrum may be expressed 



as 



G(k — ko) = A(k — ko) ■ exp 



iO(^-^o) 



aW(^ = Xo)+A» 
AW(^-^o)-exp \i^(k-ko) 



as shown in Fig. 5(b). 

Hence, substituting Eqs. (18) and (10) into Eq. (17), yields: 





-ko) 


■exp 




-ko) 




-ko)- 


exp 


i^s(k- 


-ko) 



Finally, the means of A'^^^ can be normalized to A^^\ to equalize their spectral dynamic ranges. 
Preferably, in order to avoid any degradation of the image spectrum by direct subtraction, 
o^^'' (ri) can alternatively be calculated by: 

KQ 



ifk = ko 
ifk^ko 



(18) 



(19) 



>?(n)} = 



A(k — ko) • exp 




-ko) 




A^'^(^-^o)-exp 




-ko) 


+A»5(X-^o) 



(20) 

In summary, the scattered field o^^^^ can be obtained for any illumination by Eq. (15) if an 
experimental reference field is provided. Alternatively, under the assumption of plane wave 
illumination, it can be calculated by Eqs. (19) or (20). The latter is used for processing in 
section 3. 

We compare this result to the Fourier diffraction theorem [5] of scattering potential F{fi)\ 



^{^4(^1)} = =^^)}exp [Tik,z^] . 

As) 



(21) 



It states that the scattered field uy\ recorded at plane z^, is filtered by an ideal Ewald half 

ko 



sphere = 'sj {rimk)'^ — k^ — k^ (rim- refractive index of mounting medium), and propagated by 
the latest term as known by the filtered back propagation algorithm of conventional diffraction 
tomography [6]. In our case, dividing by the 3D CTF, the spectrum is inversely filtered by an 
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'experimental' Ewald sphere. Moreover, the field propagation is intrinsically included through 
the z-dependent pre-factor in the reconstruction of Eqs. (9) and (10). Therefore, the product of 
the right side in Eq. (21) may be approximated as 

^{\(n)]='^^{o'il{n)]. (22) 

The main difference consists in the filter function. A multiplicative filter in Eq. (21) cannot cor- 
rect for realistic imaging conditions, such as aberrations or MO diffraction. It merely acts as a 
low pass filter in the shape of an Ewald sphere. By contrast, Eq. (22) implies the object function 

is) 

, which is compensated for realistic imaging conditions. In order to achieve that correction, 

fields must be divided [8]. A priori, the experimental CTF is more apt for this division since 
it intrinsically corrects for MO diffraction, apodization, aberrations and non-design imaging 
conditions, as discussed in section 3.1. 

The function F(K) is the 3D Fourier transform of the scattering potential F{r\) derived by 
the inhomogeneous Helmholz equation of the medium n{r\), where: 

n{n)=[nl-FLB{h)/k^]''\ (23) 

and n{r\) is the complex refractive index. The real part of Eq. (23) is associated with refraction 
while its imaginary part is related to absorption. 

If one were to measure the scattered field in the far zone for all possible directions of inci- 
dence and all possible directions of scattering one could determine all the Fourier components 
F{K) of the scattering potential within the full Ewald limiting sphere of 2^ = 47r/A. One could 
synthesize all Fourier components to obtain the approximation 



FLB{n) = j^ /// /^(/t)exp 



dKjcdKydK^, (24) 



{2n)^ JJjk<2k 

called the low-pass filtered approximation of the scattering potential [5] that gives rise to 
diffraction tomography. Opposed to this full approach, the approximation of the partial low- 
pass filtered scattering potential F in only one direction of gives rise to optical sectioning, as 
discussed in section 3.3. 



3. Applications 

In this section, general imaging aspects of the proposed method and their impact on phase signal 
is evaluated. Moreover, the extraction of a scattered object field is practically demonstrated to 
result in optical sectioning. 

3.1. Coherent imaging inversion 

First, the general impact of complex deconvolution on the coherent imaging inversion under 
non-design conditions is discussed. For this purpose, experimental images of non-absorbing 
mono-dispersed polystyrene microspheres {jisph = 1-59, 0 ^ 5.8/im) in water {nm,H20 = 1-33) 
are recorded at a tilt angle of 6 = 15°, shown in Fig 6(a) for raw data and in Fig. 6(b) after 
deconvolution. 

If complex deconvolution is successful, then a number of improvements in the complex field 
should be noted, accordingly indicated in Fig. 6 by regions of interests (ROI): 

Background extinction: The transparent sample images are recorded with the incident light 
o^'^ in direction of ko. According to Eq. (20), a DC value is added to the APSE to compen- 
sate for (9*^^) well seen by reduced background haze in the amplitude (cf. Fig. 6 ROI-1). 
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Fig. 6. Complex fields of polystyrene microspheres in water at a tilt angle of 0 = 15°, 
according to Fig. 2. The main images show the raw amplitude (a) and the deconvolved 
amplitude (b) with background ROI-1 and object ROI-2 with circle 0 ^ S.^fim. The inserts 
in ROI-3 show the phase parts, respectively. Colorbar, Scalebar: 4/im. 



Similarly, the removed background results in full 27r-dynamic range of phase oscillation 
(cf. Fig. 6 ROI-3). Finally, the object ROI-2 in Fig. 6 displays improved contrast since 
the objects' edges are sharpened by the complex deconvolution. 

Diffraction pattern suppression: A second motivation consists in correcting the diffraction 
pattern of the MO's APSE This correction is in particular required for high-NA imaging 
systems since the APSF diffraction pattern may result in incorrect tomographic recon- 
struction in the near resolution limit range. The diffraction pattern can be observed to be 
well suppressed by comparing ROI-1 in Fig. 6. As a result, the diffraction pattern of the 
refractive index mismatched sphere [24,25] becomes apparent in ROI-2 of Fig. 6. 

Complex aberration correction: Aberrations are intrinsic to MO, especially for high-NA 
objectives [2]. Additionally, experimental MO non-design conditions may introduce 
symmetric aberrations [3]. Those conditions include non-design refraction index of rii, 
mismatch of rii and rim, defocus of object in rim, and non-design coverslip thickness or 
refraction index. Also, axially asymmetric aberrations are introduced by the sample ro- 
tation [4, 20] . This aberration can be recognized as the asymmetric deformation of the 
diffraction pattern in the raw images (cf. Fig. 6 ROI-1 and ROI-3). As a consequence, the 
raw object in Fig. 6 ROI-2 is deformed, too. However, after deconvolution in ROI-1 of 
Fig. 6 the asymmetric diffraction pattern is removed. In the same manner, the phase oscil- 
lation becomes equally spaced after deconvolution as shown in Fig. 6 ROI-3. Eventually, 
the object can be reconstructed in a deformation-free manner in Fig. 6 ROI-2. 

Note that even with an accurate APSF and an effective deconvolution algorithm, deconvolved 
images can suffer from a variety of defects and artifacts described in detail in reference [15]. 

3.2. The impact of phase deconvolution 

The source of most artifacts is due to noise amplification as mentioned in section 2.1. The 
suggested 3D deconvolution of complex fields has the capacity to tune between complex and 
phase deconvolution according Eq. (8). Thus, noise amplification can be excluded for T = 1 
while the phase part still leads to image correction according to the previous section. As seen 
by Eq. (19), the phase deconvolution effectively acts as a subtraction of the diffraction pattern 
in phase. Strictly speaking, the recorded phase is not the phase difference between object and 
reference beam, but also includes the MO's diffraction due to frequency cutoff. The coherent 
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system is seen to exhibit rather pronounced coherent imaging edges, known as 'ringing' [26]. 
For biological samples [27], the diffraction influence in phase may be of great importance for 
the interpretation of phase signal. 

Consequently, we investigate the complex deconvolution's influence on biological samples' 
phase signal. The samples are human red blood cells (RBC) that have been fixed by ethanol 95% 
and are suspended in a physiological solution (Um^sol = 1-334 at A = 682nm). This preparation 
method allows fixing the RBC directly on the coverslip surface to avoid defocus aberrations, 
and a space invariant APSF may be assumed within the field of view [28]. The experimental 
APSF used for processing was acquired under identical object imaging conditions, i.e. same tilt 
positions with identical covers lips as shown in Fig. 2. The experimental APSF can therefore be 
accurately accessed. A comparison between the raw phase images and the phase deconvolved 
images is shown for two RBCs in Figs. 7(a) and 7(c). 



(a) RAW PROCESSED (c) RAW PROCESSED (d) T=10" 




Fig. 7. Human RBCs in phase. Images (a) and (c) show the phase images of two RBCs 
at 6 = 0. Unprocessed images are labeled 'RAW and the label 'PROCESSED' indicates 
the deconvolved phase for t = 1. The profiles in (b) and (d) compare the according height 
differences of images above, at central sections (indicated by flashes). The error bars in- 
dicate the level of phase noise O.lrad). Images (d)-(f) show the top view on RBC (c), 
processed with different T (expressed in units of Ay). Colorbars, Scalebars: 4jim. 



(a) e=2° (b) 9=4° (c) 9=6° (d) 9=8° (e) 9=10° (f) 9=14° 




Fig. 8. Inclined human RBCs in phase at various tilt angles 6. Images (a)-(f) are unpro- 
cessed and images (g)-(l) are pure phase deconvolved T = 1. According to Fig. 2, the axis 
of rotation 'x' is indicated. Colorbar, Scalebars: 4/im. 



The influence of the phase deconvolution can be seen directly by comparing these topo- 
graphic images. Based on their shapes, RBCs are classified in different stages [27]. The raw 
image in Fig. 7(c) resembles a trophozoite stage while its phase deconvolved image reveals a 
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ring structure. Similarly, the processed RBC in Fig. 7(a) reveals a trophozoite stage. A more 
detailed comparison is given by its central height profiles in Figs. 7(b) and 7(d), calculated by 
assuming a constant refractive index of hrbc = 1-39. It shows that the deconvolved phase pro- 
file follows basically the raw trend. However, similar to coherent imaging ringing, edges are 
less prone of oscillations after phase deconvolution, in particular on the RBC's edges. Thus the 
impact of the complex deconvolution is to de-blur the phase signal. 

However, the recovered phase signal is highly sensitive to noise, as demonstrated in 
Figs. 7(d)-7(f). Decreasing value of T, implies a decreased SNR through amplitude noise am- 
plification. As a consequence, random phasers can be introduced that degrade the phase signal. 
Figures 7(e) and 7(f) show that signal degradation affects most prominently RBC's central and 
border regions. For a pure retrieval of phase that is unaffected by noise amplification, T should 
not be smaller than -IdB. 

In Fig. 8, measurements under non-design imaging conditions are presented. The RBC's raw 
phase measurements become strongly deformed with increasing coverslip tilt [cf. Figs. 8(a)- 
8(f)]. This observation is in accordance with section 2.3 where strong coma-like aberrations 
are observed. Consequently, APSFs acquired under same non-design conditions, i.e. same tilt, 
can be used to effectively correct for coma aberration within the depth of field. Accordingly 
processed RBC measurements successfully recover the ring stage, even at steeper tilting angles, 
as demonstrated in Figs. 8(g)-8(l). Note that the phase signal is increased for steeper 6 since 
projection surface along y-axis decreases, as well observed in Figs. 8(g)-8(l). 

3.3. Scattered object field retrieval and optical sectioning 

From the intensity deconvolution point of view, pseudo 3D microscopy can be achieved by re- 
duction of out-of-focus haze [14], which means that the spread of objects in the z-direction is re- 
duced. Recently, this potential has also been demonstrated for digitally refocused 3D fields [13] 
and compressive holography [16]. Optical sectioning effects are therefore intrinsic to 3D com- 
plex deconvolution if T << 1 and amplitude information is not discarded. From the inverse 
filtering point of view, the z-confinement arises from the filtering by the 3D CTF as shown 
in Fig. 4. To demonstrate that scattered object field retrieval features this effect, complex de- 
convolution was performed on the RBC sample of the previous section 3.2 with T determined 
by CTF's noise level. For optimal amplitude contrast, the value of T was calculated by a his- 
togram based method, Otsu's rule, well known in intensity deconvolution [29]. The resulting 
T — 3dB, lies beneath the pure phase deconvolution criterion of -IdB, but it allows effective 
diffraction pattern suppression. Therefore, high amplitude contrast compromises phase signal. 
The results are shown in Fig. 9. 

The raw 3D field \U\'m Fig. 9(a) shows the object spread along the axial direction according 
to Eq. (9). As expected, the reconstruction features no axial confinement and the RBCs cannot 
be recognized. 

On the other hand, the scattered object field \o^^^ \ after truncated inverse filtering according 
to Eq. (20) is depicted in Fig. 9(b). It shows that the background field and out-of-focus haze 
are successfully removed. The RBCs' edges can be identified as strong scattering objects and 
the scattered field can be recognized to match in size and position the anticipated RBC values. 
Although the image quality is affected by artifacts in axial elongation [30], its axial dimen- 
sions match well. However, the refractive index mismatched cell membrane {nupid > 1.4) and 
mounting medium (nm^sol = 1-334 at A = 682^m) result in strong scattering well visible in the 
xz sections. 

Finally, the fields related to refraction, \n\, can be reconstructed in Fig. 9(c) according to 
Eq. (23). In particular, the refraction due to the strong scattered field around the RBCs allows 
a good three-dimensional localization of the RBCs' edges. Moreover, the higher refraction 
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Fig. 9. 3D rendered images of two human RBC at 0 = 0. Images (a), shows \U\ in 3D- 
space in the middle row. Bottom and top rows show the sections through central RBC 
positions (indicated by flashes). Accordingly, the field lo^^^^ \ is represented in (b), and \n\ in 
(c) (uncalibrated levels), respectively. The RBCs' positions are compared to an oval area of 
6fim X 2.5 fim, based on measurements from section 3.2. Colorbars. 



index due to its hemoglobin content is well visible in the xz sections. Note that these data 
are reconstructed from only one hologram for a single incident angle and therefore missing 
angles in Eq. (24) affect the reconstruction as seen by the lateral artifacts. However, if the 3D 
inverse filtering technique is combined with a multi-angle acquisition, it holds the potential of 
quantitative 3D refraction index reconstruction. 

4. Concluding remarks 

In this paper, we have described theory, experimental aspects as well as applications of a novel 
method of 3D imaging using realistic 3D CTF inverse filtering of complex fields. 

Our theory connects three-dimensional coherent image formation and diffraction theory and 
results in a model for object scattering reconstruction by inverse filtering. This approach is ex- 
perimentally complimented by the ability to characterize the DHM setup with background free 
APSF thanks to the use of a complex point source. The physical importance of the realistic 
3D CTF is demonstrated with experimental data. The technique features effective correction of 
background illumination, diffraction pattern, aberrations and non-ideal experimental imaging 
conditions. Moreover, the regularization of the three-dimensional deconvolution of complex 
fields is shown to yield reconstruction in the complex domain. Depending on the threshold, 
phase de-blurring or optical sectioning is demonstrated with RBC measurements. Most essen- 
tially, the capability of scattered field extraction is experimentally presented. 

In conclusion, the demonstrated technique bears the potential to reconstruct object scattering 
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functions under realistic high-NA imaging conditions which play a key role in high resolution 
diffraction tomography. 
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